The Cosmic Microwave Background and Helical Magnetic Fields: 

the tensor mode 



Chiara Caprini/'0 Ruth Durrcr,^'Q and Tina Kahniashvili^'0 

^Department of Astrophysics, Denys Wilkinson Building, Keble road, Oxford 0X1 3RH, UK 
^ Departement de Physique Theorique, Universite de Geneve, 
24 quai Ernest Ansermet, CH~1211 Geneve 4, Switzerland 
"'Department of Physics and Astronomy, Rutgers NJ State University, 
136, FreUnghuysen RD., Piscataway, NJ, 08854-8019, USA 

and 

Center for Plasma Astrophysics, Abastumani Astrophysical Observatory, 2A, Kazbegi ave., Tbilisi, 380060, Georgia 

(Dated: February 2, 2008) 

We study the effect of a possible helicity component of a primordial magnetic field on the tensor 
part of the cosmic microwave background temperature anisotropies and polarization. We give an- 
alytical approximations for the tensor contributions induced by helicity, discussing their amplitude 
and spectral index in dependence of the power spectrum of the primordial magnetic field. We find 
that an helical magnetic field creates a parity odd component of gravity waves inducing parity odd 
, polarization signals. However, only if the magnetic field is close to scale invariant and if its heli- 

' cal part is close to maximal, the effect is sufficiently large to be observable. We also discuss the 

' implications of causality on the magnetic field spectrum. 
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I. INTRODUCTION 



The observed Universe is permeated with large scale coherent magnetic fields. It is still under debate whether these 
magnetic fields have been created by charge separation processes in the late Universe, or whether primordial seed 
1^ ' fields are needed. Recently, it has been proposed Q that also 'helical' magnetic fields, i.e. fields with a non- vanishing 

■ component in the direction of the current, B • (V x B) 7^ 0, could be produced e.g. during the electroweak phase 
' transition (see also 0). 

Extended studies have already investigated effects of stochastic magnetic fields with vanishing helicity on the cosmic 
microwave background (CMB) (see 3, 4,^,^] and others). In a seminal paper •?,], Pogosian and collaborators have 
\ investigated the possibility that a helical magnetic field can induce correlations between the temperature anisotropy 

■ and the B mode CMB polarization. 

Qh[ In this paper we want to go beyond that work. We determine all the effects on the CMB induced by a helical 
Q magnetic field. We shall actually show that, contrary to the statement in Ref. 0, a helical component also introduces 
%-[ ' pure CMB anisotropies and polarization. But of course its most remarkable effect is the above mentioned correlation 
CZ3 . of temperature anisotropy and B polarization. We shall show that also a correlation between E and B polarization 
' is induced. 

^ , In this paper we discuss only the tensor mode, gravitational waves, since the calculations for this case are simplest. 

Even if the resulting observational effects are small and may not be detectable, we find it interesting since it is 
^ . completely new and contains several surprising elements. Furthermore, a fluid vorticity fleld or non parity invariant 
initial spectrum of gravitational waves produced during inflation could induce very similar effects; in that sense our 
results are more generic than their derivation. 

In the next section, we discuss the magnetic field spectrum and define its symmetric and helical contributions. 
Then we compute the tensor component of the magnetic field energy momentum tensor which acts as a source for 
gravity waves. In Section IV we determine the induced gravity wave spectrum which also has a symmetric and a 
helical contribution. In Section V we compute the induced CMB temperature anisotropy and polarization spectra as 
well as the above mentioned correlations. Finally, we discuss our results and draw some conclusions. The paper is 
complemented by an appendix where details of calculations and tests of some approximations can be found. 
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II. THE MAGNETIC FIELD SPECTRUM 



We consider a primordial stochastic magnetic field created before equality, during the radiation-dominated epoch (or 
earlier). During this period of the evolution of the Universe, the conductivity of the primordial plasma on scales larger 
than the Silk scale A > As is very high, effectively infinite 3- Hence, the 'frozen-in' condition holds, E = — v x B, 
where v is the plasma flux velocity, E is the electric field induced by plasma motions and B is the magnetic field. 
Moreover, large scale magnetic fields always induce anisotropic stresses, so that their energy density i?^/87r must 
be a small perturbation, in order not to break the isotropy of the Friedmann Robertson Walker background. This 
allows us to apply linear perturbation theory. Both, the magnetic field energy and the plasma peculiar velocity are 
treated as first order perturbations; consequently, the energy density of the induced electric field will be 3rd order in 
perturbations theory, and can be neglected. Also terms EiBj are of second order and therefore neglected. 

At sufficiently large scales, it is possible to neglect the effects of back reaction of the fluid on the evolution of the 
magnetic field: the time dependence decouples from the spatial structure, and, due to flux conservation, the magnetic 
field evolves like B(77, x) = B(?7o, x.)/a{r])'^, where we use the normalization 0(770) = 1 and a subscript denotes today. 
At smaller scales however, the interaction between the fluid and the magnetic field becomes important, leading mainly 
to two effects: on intermediated scale, the plasma undergoes Alfven oscillations, and B^{k) B^{k) cos^ {vAkrj) (where 
v\ = B^/{4:Tr{p + p)) is the Alfven velocity, here B is the field averaged over a scale of order vat]); on very small 
scales, the field is exponentially damped due to shear viscosity 0,0,0i^3- -R-^f. 0, we will account for this 

damping by introducing an ultraviolet cutoff at wavenumber k£){rf) in the spectrum of B (see also 6]). 

Following Refs. jlj,t3], we introduce an hclicity component A{k) in the magnetic field two point correlation function: 

(B,(k)Br(k')) - ^^'5(k - \<^)[PjiS{k) + ie,i^kmA{k)] , (1) 

where S{k) and A{k) are respectively the symmetric and helical part of the magnetic field power spectrum. Pij = 
5ij — kikj is the usual transverse plane projector satisfying the conditions PijPjk — Pik, Pijkj = 0, e^/ is the totally 
antisymmetric tensor, and ki = ki/k. We use the Fourier transformation convention 



Bj{k) = J d^x exp(ik • x)Bj(x), Sj (x) = -j^^ j dPk exp(-ik • x)Bj (k) 



(2) 



For simplicity, as in Refs. 0, Q and others, we shall assume that the magnetic field is a Gaussian random field. 
Then all the statistical information is contained in the two-point correlation function and the higher moments can be 
obtained via Wick's theorem. 

As explained in Ref. 01, the magnetic field helicity is determined by (B-(VxB)). For a better physical understanding 
of the effects which this new helicity term has on CMB anisotropics, it is useful to introduce the orthonormal 'helicity 



ty ter] 

basis' (e+, e~, 63 = k) (see also 0,^2), where 



e±(k)---^(ei±»e2) , (3) 

and (ei, 62, = k) form a right-handed orthonormal basis with 62 = k x ei. Under the transformation k — k 
we choose 62 to change sign while ei remains invariant. The basis (e+, e^, k) has the following properties: • = 
— 1, • = 0, and e='=(k) = e^(— k), as well as ik x e='= = ±e^. The components of a vector with respect to 
this basis will be indicated by a superscript ±. For a fixed (fc-independent) basis we will instead use the usual Latin 
letters as indices. An arbitrary transverse vector v can be decomposed as v = i;+e+ -|- v^e^ . Here w"*" is the positive 
helicity component and is the negative helicity component. 

With the definition and the reality condition (i3^(k))* — — i3*(— k), we obtain the connection between the 
power spectra 5(/e), A{k) and the magnetic field components in the new basis: 

- (S+(k)B+(-k') + B"(k)B-(-k')) = (27r)^5(fc)5(k-k') , (4) 
(S+(k)B+(-k') -B"(k)B-(-k')) = (27r)3yl(fc),5(k - k') . (5) 

In other words, A{k) represents the difference of the expectation values of the positive and negative helicity field 
components. If A does not vanish, the left handed and right handed magnetic fields have different strength. 

We assume that both the symmetric and helical terms of the magnetic field power spectrum can be approximated 
by a simple power law 0: 

Sik) = { ^'^ k < ko /gN 

^ ^ 10 otherwise ^ ' 



and 

A(k) = Mo^"^' forfc<fcz5 
^ ' 1 otherwise ^ ' 

where Sq, Aq are the normahzation constants, and ns, ua the spectral indices of the symmetric and hehcal parts 
respectively. 

With (0 d) , we can express the normalization constants 5*0 and Aq in terms of the averaged magnetic field energy 
density Bx^ = (B(x) • B(x))|a, and the absolute value of the averag ed helicity = A|(B(x) • (V x B(x)))||a 
respectively, both smoothed over a sphere of comoving radius A. B\ measures the amplitude of helicity on the given 
comoving scale A. 

In order to calculate these quantities, we convolve the magnetic field and its helicity with a 3D-Gaussian filter 
function, so that Bi Bi * />, where f\{k) — exp(— A^fc^/2). The mean-square values B\ and B\ are then given 
by the Fourier transform of the products of the corresponding spectra S{k) and kA{k) with the square of the filter 
function fy. 



?2 _ ^ /■ ./.M ? /,x2 l^ol 1 ^fnA 



In order not to over-produce long range magnetic fields or helicity as /c — > 0, we require for the spectral indices 
ns > —3 and ua > —4 (for ua < —3 and ua < —4 the integrals ^ and lO diverge at small k). 

Using © , (13 and the definition of the magnetic field spectrum , we can rewrite expressions (0J and © in the 
form (see also j3) 

-{B+ik)B+i-k') + B-ik)B-i-k')) = i2nf ^^l%^ {XkrSik-k'), (10) 

{B+{k)B+{-k')-B-{k)B-{-k')) i2nf-^-^{Xkr-S{k-k') , (11) 

for k < kjj and for k > ki). 
Using that 



hm |((k X B(k)) . B(-k'))| < Jim (B(k) • B(-k')) 

k'^k k'^k 



we can conclude that 



S{k) > \A{k)\ . (12) 

Since S{k) (x (|Bp), it is clear that S{k) > 0. The reality condition requires Aq to be real, but it can be either 
positive or negative. For Eq. H12I) to be valid on very small values of k requires 

nA>ns. (13) 

Applying Eq. l(T^ also close to the upper cutoff k^, we have in addition 

\Ao\ < S^kl'-^- . (14) 

In terms of the magnetic fields on scale A this gives roughly 

Bf < BlikoXT'-''^ . (15) 

Usually the damping scale is much smaller than the physical scale of interest, A so that Xk^ ^ 1. Therefore, if 
ns — nA ^ 0, the helical contribution is significantly suppressed on all scales A > A^ = 1/ku- As wc now show, this 
is always the case if the magnetic field is causally produced. 

Most mechanisms to produce magnetic fields with a helical component are causal. By this we mean that all 
correlations above a certain scale, usually some fraction of the Hubble scale at formation, have to vanish. If this 
is the case, causality implies an additional interesting constraint, which we now derive. For this we assume that 
the correlation functions {Bi{x)Bj{y)) and (_Bi(x)(V x B(y))j) have to vanish for |x — y| > R for some scale R. 



Hence they are functions with compact support, which imphes that their Fourier transforms, PijS{k) and eijikiA{k) 
are analytic functions. Therefore, for sufficiently small values of k they can be approximated by power laws as in 
Eqs. (j6l7|l . Since fcj is not analytic but kkj is, this implies 

ns>2 and > 1 , (16) 

where ns has to be an even integer while ua has to be an odd integer. But since we need ua > ns^ this leaves us 
with 

'"^s ^ 2 , an even integer and (17) 
UA > 3 , an odd integer. (18) 

Causality together with the condition H12(l leads to an additional suppression of helical fields on large scales. Also 
ordinary causal magnetic fields cannot be white noise but are severely suppressed on large scale due to the non-analytic 
pre-factor Pij in the power spectrum which is a simple consequence of the fact that magnetic fields are divergence free 
V • B = 0. This has already been discussed in Refs. 0,101 . The causality constraint need not to be satisfied if the 
magnetic fields are generated before or during a period of inflation where the causal horizon diverges. For a detailed 
discussion of causality see [T3|. 



III. MAGNETIC SOURCE TERM FOR TENSOR METRIC PERTURBATIONS 

The anisotropic stresses which act as source for metric perturbations are given by the magnetic field stress tensor 

(k) = ^ ^ j Sp mv)B* (p - k) - (p)i3r (p - k)%] . (19) 

Here we are interested in the generation of gravitational waves, and consequently we need to extract the transverse 
and traceless part of Tij. The form of a general projection to extract any mode (scalar, vector or tensor) from a 
generic tensorial perturbation can be found in [l5l |. We make use of the tensor projector Tijim — PuPjm — \PijPim 
(see also 0). The tensor contribution to Tim is given by 



Hy = {PaPjm - T;PvPlrn)Tlm ■ (20) 



Moreover, since the magnetic field is a stochastic variable, we need to calculate the two point correlation tensor of 
Ty(k), which takes the form 

(T,,(k)rr„(k')) = J^J^ I d'p I d\ (S.(p)S,(k - p)BK-q)S™(q - k')) +•••% + ••• 5im , (21) 

and we are not interested in terms proportional to Sij and dim , which after being projected out will not contribute to 
the final result for the tensor perturbation (HyH/m) (see appendix A. in Q). Before applying the tensor projection, 
we can simplify the right hand side of H21|l using Wick's theorem, expressing the four point correlators in terms of the 
two point ones, 

{B,{k,)B,{kj)Bi{ki)Bm(km)) = {B,ik,)Bj{k,)){Bi{ki)Bm(km)) 

+ {B,ik,)Bi{ki)){B,{kj)Bmikm)) 

+ {B,{k,)Bm{km)){Bj{kj)Bi(ki)) . (22) 

Since the two point correlation function given in Eq. |^ is not symmetric, we are not allowed to change the order 
of indices i,j,l,m inside an expectation value. With Eq. JQ) we can then compute the correlation function H21|) 
which consists of a purely symmetric part proportional to J d'^pS{p)S{\k — p|) , a purely helical part proportional to 
/ (PpA{p)A{\k — p|) , and mixed term, i J (PpS{p)A{\k — p|) (the full expressions are given in Appendix A, Eq. (^3). 
The first two terms contribute to the symmetric part of the two point correlation function of the tensor source, while 
the two latter terms give rise to a helical contribution. To express them we now introduce the two point correlation 
function for the tensor source, which can be parameterized as 

(H,,(k)Hr„.(k')) ^ ^ [M.,imfik) + zA^^.imgm ^(k - k') , (23) 



where the tensors Mijim and Aijim are given by 

Clearly, both Mijim and Aijim are symmetric in the first and second pair of indices. Mijim is also symmetric under 
the exchange of ij with Im while Aijim is anti-symmetric under this permutation. We shall often use simple properties 
like 

M^jiJ = 4 , MiUrn = MijH = (26) 
PqiJ^ijlm — -^^qjlm •, Pqi-^ijlm — •^qjlm (27) 
•^^ijlni-^^ijlm ^ ^ijhri'^ijlm ~" 8 (^^) 
'^ijlm-^^ijlm — i '^ij'ij ~" ^iijl " ^ijll ^ ■ (^^) 

According to Eq. (I20|) . we have now to act on (rafc(k)r^^(k')) with the tensor projector 

^"'^(k, k') = {P^aPJb - IP.jPabmiPlcPrnd - ^P;mPcd)(k') • (30) 

In these calculations we don't need to care about the position (up or down) of Latin indices as they are always 
contracted by a Kronecker S. The symmetric and antisymmetric parts of Eq. H23|l are invariant under the application 
of the projector (|30|l . so that it is easy to separate the symmetric and helical parts of the source spectrum, /(fc) and 
9ik): 

{27rfS{k-k')f{k) = ^MabcdiTabikKA^')} (31) 

(2^)35(k-k')5(fc) = ^Aabcd{rabik)T:,ik')) . (32) 

Moreover, by applying the tensor Mijim to Eq. HA1|I of Appendix A, we obtain (the first term of this has already 
been computed in Refs. 0,|E0) 



y d^p[5(p)5(|k-p|)(l+7^)(l + /3^) + 4A(p)A(|k-p|)(7/3)] , (33) 

where 7 = k • p and /3 = k • (k — p). Note that the square of the helical part of the magnetic field spectrum 
contributes to the symmetric part of the source spectrum. This is not surprising, since the product of two quantities 
with odd parity has even parity. The antisymmetric part of the source spectrum is obtained by acting with Aijim on 
Eq. IjAip of Appendix A. It is given by the mixed terms, 

/ d'p Sip) A{\k - p\){l + . (34) 
We can also express the correlator (|23ll in terms of the basis efj introduced in [lll |. 

efj = - Y ± ie2)t x (ei ± 162)^- . (35) 

These form a basis of tensor perturbations, satisfying the transverse-traceless condition Sijefj = 0, kicfj = and 
efj^fj = 3/2. Positive circularly polarized gravity waves are proportional to e^, while negative circularly polarized 
gravity waves are given by the coefhcient of e~- . In this basis Ily is expressed as 

%(k) = ejn+(k)+e-.n-(k) . (36) 

We can rewrite f{k) and g{k) in terms of the components 11^ as 

(27r)3j(k - k') f{k) ^ (2^)3^(k - k') |n(fc)|2 = ^ (n+(k)n+*(k') + n-(k)n-*(k')) , (37) 

(27r)3j(k-k') g{k) = (n+(k)n+*(k')-n-(k)n-*(k')) . (38) 



Here we have used the form of A4 and A in this basis, 



•Aijlni 

and the simple properties of Mijim and Aijim mentioned above. One also has 



3 L 

4i 

J [4j ® e,-;„ - e-- (g) e+„] , 



(n+(k)n-(k')> 



(27r) 



<5(k-k') (/(fc)+5(fc)) 



Similarly, defining the usual linear polarization basis 



= (ei X ei - 62 X 62) 
= (ei X e2 +62 X ei)^^. 

and the components of 11 with respect to this basis, 
we obtain also 



(n'^(k)n'^*(k') + nx(k)nx*(k')> = {2nf 6 {k-k') f{k) 
(nx(k)n^*(k') -n'^(k)nx*(k')) = i{27rfs{k-k')g{k) 



With Eqs. 1221 El, we find 

(27r)3[/(fc)+g(fc)] = 



(47r); 



1 d^p [S{p){l + 7^) + 2A{ph] ■ [Si\k - p|)(l + P^) + 2A{\k - p\ 



Let us introduce the tensor 



so that 



(27r)3 47r 



(47r) 



[Py(k)5(fc)+zey,fc,A(fc)] 



;i?,(k)B*(k')) =^(k-k')Q.j(k) ; 



with Qij (— k) = Q*^ (k) one then finds 

f{k)+g{k)^ P.ij{k) ~ ieijqkq Pim{k) + if-irnq'kq' / d^p (5y(p)(3^„,(k- p) 



(39) 



(40) 
(41) 



(42) 
(43) 



(44) 



(45) 



(46) 



(47) 



Using Eqs. 1)619(1 . 1)33(1 and 1(34(1 . it is possible to calculate f{k) and g{k). The details of the calculations are given 
in the Appendix A. The integrals cannot be computed analytically, but a good approximation gives, for k < ko (see 
also nil): 



f{k) ~ As {Xko) 



\2ns+3 



ns 



ns 



-(Afc) 



2ns+3 



g{k) :^ aXkn) 



(Afc) 



nA - 1 
ns + i \k 



- A A (Xkn) 



ns+nA+2 



nA + 4 



(Afc) 



(48) 
(49) 



where As, Aa and C are positive constants given in Eqs. l(A13p to ((Af 5p of Appendix A. They depend on the 
spectral indices ns and nA of the magnetic field and on its amplitudes, which are given in terms of _B|, Bj^, and A. 

Note that the contribution of magnetic field helicity to the symmetric part of the source, /(fc), is negative. But 
it is easy to check that Eq. 1(12(1 insures that it never dominates, hence / > 0. For ns,nA > —3/2, the two terms 
proportional to the upper cutoff k'^^'^'^^^ dominate in /(fc), which consequently depends only on the cutoff frequency 
and behaves like a white noise source (4|. For ns < —3/2 or also nA < —3/2, the dominating terms go like 



and respectively. On the contrary, the antisymmetric source g{k) never shows a white noise behavior. For 

ns + nA > —2 the dominant term is proportional to fc fc^f For ns + ua < —2, g{k) does not depend on the 

upper cutoff, but is proportional to fc"s+nA+3 ^j-^g singularities in the pre- factors As, A a and C which appear at 
ns = —3 and ua = —4 are the usual logarithmic singularities of scale invariant spectra. But as mentioned in Section II 
the helical contribution must obey ua > ns > —3. The apparent singularities in the pre-factors at ns^A = —3/2 and 
at ns + nA — —2 are removable when multiplied with the fc-dependcnt parts as in Eqs. H48|l and 149|) . In the integrals 
over k which we shall perform to calculate the CiS we only take into account the dominant terms. 
If the magnetic field is causal, we expect ns = 2 and tia — 3, so that 

/(fc) ~ AsikoXy - AAikoXf (50) 
g{k) ~ CkXikoXy . (51) 

Comparing the limit given in Eq. (O with the expressions for As and Aa derived in the Appendix A, it is easy to 
see that / always remains positive. 

The analysis of the evolution of a non-helical magnetic field interactiiig with the primordial plasma, and the 
derivation of the appropriate damping scale kjj, has been discussed in Refs. Q and [lfl| |. where the authors considered 
a magnetic field with a tangled component superimposed on a homogeneous field. We assume that the latter can be 
obtained by smoothing our stochastic field on a scale which is larger than the damping scale (for details, see |^[T3|). 
The damping scale for the tensor mode is obtained taking into account that the source of gravitational radiation 
after equality becomes sub-dominant so that the relevant tensor damping scale is the Alfven wave damping scale from 
the time of the creation of the magnetic field up to equality |l2j . Since we are interested here in the imprint of the 
magnetic field on the CMB, we need not to care about the time evolution of the damping scale, the relevant scales for 
the CMB tensor anisotropics being those which are greater or equal to the horizon at equality. Therefore, the relevant 
cutoff scale is given by the Alfven wave damping scale at equality A;^^ ~ VAl-yiTcq), where l^{Tcq) « 0.35 Mpc is the 
comoving diffusion length of photons at equality (here we have used that 1^^^^{T) ~ 10^^cm(T/Tdoc)~^, from as 
well as Zcq — 3454 and z^cc — 1088 from the WMAP results 0|). The Alfven speed is at most of order 10"'^, so that 
the damping scale is on the order of kpc or smaller. 

Even if considering an helical component in the magnetic field, we set all the power to zero on scales smaller than 
kjj . This is not really correct since simulations show 

m that the spectrum dimply decays like a power law with 
index of the order of —4 on small scales, k > k^. However, as we shall see, for ns.A < —3/2 the induced C^'s are 
dominated by the contribution at the largest scales, fc^j^, for the kinks, ns^A ^ —4 part of the spectrum. Therefore, 
we do not loose much by neglecting the contribution from the scales smaller than fc^^. 



IV. MAGNETIC FIELD INDUCED TENSOR METRIC PERTURBATIONS 



A stochastic magnetic field can act as a source for Einstein's equations and hence generate gravitational waves, see 
for example [3. 1^ II^. The tensor modes are the simplest case of metric perturbations, and in the transverse and 
traceless gauge they are fully described by the tensor /iij(x, 77), satisfying 

= hji, ha = 0, hi^y = . (52) 
The linear evolution equation for gravitational waves is 

/iy(k,?7) + 2-Ay(k,7?) + fc2/iy(k,77) = 4^n,j(k), (53) 

where nij(k) is the source tensor given in H20() . and we have multiplied in the time dependence 0^^(77), which comes 
from the fact that the magnetic field is frozen in the plasma. Therefore, ny(k, ry) is a coherent source, in the sense 
that each mode undergoes the same time evolution [l3| • We neglect other possible anisotropic stresses of the plasma 
(collisionless hot dark matter particles or massless neutrinos have anisotropic stresses which do source gravitational 
waves, but this effect is very small 0). 

We want to compute the induced CMB anisotropics and polarization (see Section V), which can be expressed in 
terms of the two-point correlation spectrum (/i^ (k)/i;„j(k')), taking the form jl.ll^: 

(/iy (k,r7)/iL(k',77)) - ^[A^y,„,i7(fc,77)+a.,,™H(A:,77)]<5(k-k') . (54) 

Here {2^:YH{k.r|)5{\^— k') = (/iij(k)ft,*^ (k')) is the usual isotropic part of the gravitational wave spectrum which is 
sourced by /(fc), and 7i{k,r]) describes the helical part, sourced by g{k). 



The perturbation tensor hij can also be expressed in terms of the basis ef, defined in Eq. H35(l : 



hijiK v) = h+{k, 77)e+- + h (k, 7j)e^j . 
Just like for the anisotropic stress power spectra, we now find that 

(2^)35(k-k') i?(fc,r?) EE ^(A+(k,r;)/i+*(k',7?) + /i-(k,7?)A-*(k',,?)) , 
(27r)35(k-k')7i(fc,r7) ee -^{h+{k,7^)h+* {k' ,7j) ~ (k,r^)h-* (k' ,7^)) . 

In terms of /i'^ and , defined like in Eq. (|40|l . H parameterizes the correlation between and 

(/i^ (k)/i^*(k') - /i'^(k)/i^*(k')) = i(27r)35(k - k') . 
The evolution equation for the components (k, 77) is simply 



/i±(k,77) + 2-A±(k,r;) + fc^/i±(k,?7) 



SttG 



n±(k) 



(55) 

(56) 
(57) 

(58) 

(59) 



We need to determine the functions /i^(k, 77) (see Eq. (|66f) below). An approximate solution to the above differential 
equation can be found in j3| or The important point is that because of the rapid falloff of the magnetic field 

source in the matter dominated era, perturbations created after equality (jyoq) are sub-dominant, so that one obtains, 
for the dominant contribution at 77 > 77oq: 



eq 



n±(k) 



j2ikri) 



(60) 



where fir is the radiation density parameter today and ^in.oq correspond to the redshifts at the moment of creation of 
the magnetic field and at matter radiation equality respectively. The function j2 is the spherical Bessel function . 
The term ln(zin/zcq) accounts for the logarithmic build up of gravity waves from z-^ to Zcq. For the spectra (j56|l and 
(|57jl we then obtain 



H{k,7^) 





■ 167rG / 










VZoqV 


V 




■ 167rG , / 










V^oqV 


V 



1 2 



fik) , 



(61) 
(62) 



The gravity wave power spectra H / pr and 7i/ Pr are constant on large scales, krj <^ 1 and decay and oscillate inside 
the horizon. 

Our first result is that a helical magnetic field induced a parity odd gravity wave component. From Eq. H59(l it is 
clear, that such a component is introduced whenever there are parity odd anisotropic stresses. It could in principle also 
be detected directly, via gravity wave background detections experiments. We do not discuss this very hypothetical 
idea any further, but calculate the effect of such a component on CMB anisotropics and polarization. 



V. CMB FLUCTUATIONS 



Magnetic fields in the universe lead to all types of metric perturbations (scalar, vector and tensor, for more details 
see [^). In it is shown that vector and tensor perturbations from magnetic fields induce CMB anisotropics of 
the same order of magnitude. In this paper we estimate CMB fluctuations due to gravitational waves induced by a 
stochastic magnetic field, the spectrum of which contains an helicity component, ^(fc) ^ 0. Since the CMB signature 
of chaotic magnetic fields with only an isotropic spectrum is given in detail in Refs. ^ila, here we concentrate on the 
effects from the helical part of the magnetic field spectrum, and we will discuss the corrections which it induces to 
the previous results. 

To compute the CMB fiuctuation power spectra we use the total angular momentum method introduced by Hu 
and White jTll |. By combining intrinsic angular structure with the spatial dependence of plane- waves, Hu and White 



obtained integral solutions for all kind of perturbations. The angular power spectrum of CMB fluctuations can then 
be expressed as 



c, =- dkk Y: ^T^n 27T^' (^^^ 

m— — 2 

where X takes the values of Q, temperature fluctuation, E, polarization with positive parity, and B, polarization with 
negative parity, for each perturbation mode. The index m indicates the spin, and for tensor modes m = ±2. Since 
we only consider tensor modes in this paper, we suppress the index 2 and just denote the two states by + and — in 
what follows. The description given in Ref. [a| applies the total angular momentum method to parity even magnetic 
field spectra: in this case, according to parity conservation the sum over ± can be replaced by a factor 2. In our case 
instead, we always need to sum over both states. 

From the form of /(fc), the parity even CMB fluctuation correlators can be expressed as: 

X X' 

where C'(yi)^ is the power spectrum induced by the purely helical part of the source term, proportional to 
A(p)A(|k — p|). The contribution of this helical part to the parity even CMB power spectra is always negative, 

— . X X X X 

but, as we shall see, the condition H12|) insures that Cj.^'^^ < C'j.^j^ so that the power spectra do not become negative. 

The new effect is that the helical part of the magnetic field now also induces parity odd CMB corr elators, Cf ^ 
and Cf^ (see also 7,]). These are expressed in terms of the helical magnetic source g(\i) which is proportional to the 
convolution of A{k) with S{k) (see Eq. (E^). 

We now derive the CMB fluctuations Qf{'r]o,k), E^{riQ,k), B^{r]Q,k) and then perform the integral H63|) . Rather 
than a numerical study, we present analytical approximations for our results. These are not very accurate, but allow 
a discussion of the dependence of the correlators on ns and ua- We will also be able to determine the spectral index 
of the CMB correlators (dependence on £) as a function of ns and ua- At the present stage, we think this scaling 
information is more interesting than accurate numerical results. These can than follow for specific, interesting values 
of the spectral indices in future work. For a magnetic field with no helical component, this program has been carried 
out in Ref. 0, and we shall just refer to their results but not re-derive them here. 

Below, we shall always work in the approximation of 'instant recombination'. Moreover, in our approximations 
we didn't take into account the decay of gravity waves for modes which entered the horizon before decoupling. Our 
results therefore will be reasonable approximations (within a factor of two or so) only for ^ ^ 60, where the tensor 
CMB signal is largest. Even though, this may seem poor accuracy, here we only want to obtain estimates of the 
correct order of magnitude of this anyway small effect. This will enable use to judge for which cases a more involved 
numerical study is justified. 



A. CMB temperature anisotropies 

Within the instant recombination approximation, gravitational waves simply cause CMB photons to propagate 
along perturbed geodesies from the last scattering surface to us. The induced CMB temperature anisotropies are 
given by [lOl 

e(?7o,k,n)~ / drj ejip{-i{ijo ~ v)'^- n)hij{k,ri)nihj . (65) 
In the total angular momentum formalism this becomes 



no 

dvh^{k,rj)j^[ki7jo-ri)] , (66) 



»?do, 



2^+1 3 

where are the tensor temperature radial functions of the two different parities, both given by 



The somewhat unusual factor 4/3 comes from the fact that this formula takes into account polarization, while Eq. H65() 
does not. A detailed derivation can be found in Ref. Illl. 



Using the solution H60|l for /i*(k, 77), we obtain 
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(68) 







where we have set x = krj and xq = kr/o- For the second ~ sign we have used the approximation ljB5|l given in 
Appendix B for the integral over x. This approximation is valid only for Xdcc = kri^^cc ^ 1- 
The general expression ()63|l for the temperature anisotropy power spectrum now gives 
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A good approximation for the function f{k) is given in Appendix A, Eq. ljA9|l . The first term of l|A9|l comes entirely 
from the non- helical component Bx, and has already been determined in Refs. (0,01); the second term comes instead 
from the helical component, and its influence on the Ce is new. We denote it by C^^^. Then, splitting the induced 
temperature anisotropy power spectrum as 



(70) 



we obtain (now xq is renamed a;) 
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ee 
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where we have set xd — kjjrjQ. We have introduced the 'helicity density parameter' ^Ia defined by 
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and analogously we will use 
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(73) 



where we have introduced Bf^ — S^(fc£)A)"'^+^, the field strength at the cutoff scale l/fc_D, and correspondingly for 
Bk[, ■ With these definitions the results will be expressed entirely in terms of physical quantities and the reference 
scale A does no longer enter. 

Remember also that (27r)'*(S^A"'*+'^)^/r^ ( "V'^ ) ~ l^ol^' where \Ao\ is the normahzation of the helical component 
of the magnetic power spectrum The integral (|69|l is dominated at xo — I- With xo/xdec = Vo/Vdec — 60, this 
means that our approximation is valid for ^ < 60. 

If ua > —3/2, the first term in the square bracket in Eq. (|71|l dominates. Since the integral converges and is 
maximal around k ~ ^/??o ^ kj), wc can replace it by the integral to infinity and use Eq. I)B7|I of Appendix B. This 
gives 
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(74) 



for UA > —3/2 . 



The temperature power spectrum has the well known behavior of C^'s induced by white noise gravity waves, Ci cx £. 
If UA < —3/2, the second term in the square bracket of Eq. (|71|l dominates, and we find 
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(75) 



for 



3<nA< -3/2 



Like for the symmetric contribution given in Refs. we get a scale-invariant spectrum for ua = —3. The 

expressions for i'^C^^^ are obtained from those given above upon replacing by fls, tia by ns and ( ) 
p2 (^!i£±2.y For -3 < ns < -3/2, one also has to replace the factor (ua - + 4) by ns/{ns + 3). We do not 

repeat these formulas here since they can be found in Ref. (up to some factors of order unity which are of no 
relevance for this discussion). 

This is in principle the final result for temperature anisotropics. Let us check that C^^^ is indeed never larger than 

Cfg^^ so that 



We first consider ua > ns > —3/2. Then 

^ Bt r^(^) {2ns + 3) (fcgA)^("^-"^) ^ [App 2^5 + 3 
C©© Bl T^^) {2nA + 3) sl fc^*"^-"-^) + 3 



< 1 ■ (76) 



In the first equality we have inserted the definitions of and and the last inequality comes from Eqs. H14|l 
and H13|l. If instead ns < ua < —3/2, we find 

-iV(nA,ns)— ^7— — T , (77) 



where N(nA, ns) is a function of the spectral indices ns and n^. It is of order unity in the allowed range, —3 < nA < 
ns < —3/2. Now knTjo ^ £ for all values of £ for which our result applies. Hence again 

(J00 

7^ < 1 • (78) 

Finally, we consider the case —3<ns< —3/2 < nA, so that we have to apply the result ifT^ for C^A)e ^^^^ f?^ 
with the mentioned modifications for C^J^^. A short calculation gives 

since the first factor is less than one due to Eq. p4|l and kjj ^ £/riQ with ns < —3/2. 
Clearly, the helical component is maximal for nA — ns, where we may have |^o| — Sq. 



B. The induced CMB polarization 



Tensor perturbations induce both E polarization with positive parity, and B polarization with negative parity. 
CMB polarization induced by gravity waves has been studied for example in Refs. JJ, 21, 22J, while the contribution 
from a magnetic field has been discussed in l^^l • Our aim is to estimate the effect on the polarization signal from 
the helical component of the magnetic field. Like for the temperature anisotropics, we use the angular momentum 
method developed in Ref. [lll |. 



1. E type polarization 

The integral solution for E type polarization from gravity waves is given in |lll |. Again, we will work in the 'instant 
recombination' approximation. The order of magnitude of our result is still reasonable for £ < 60, since in this case 
also we restrict ourselves to the evaluation of the super-horizon scales spectrum. In our approximation we have 

^^TT^ = ^j\ P dr;A±(k,ry).±[M%-r7)] , (80) 



here 



4W = 7 



for ^ > 1 



(81) 



is the E-type polarization radial function for the tensor mode , and for the last equality we have used the recurrence 
relations for spherical Bessel functions (|B14I IbT5|| . 

We now use our solution ^^^i to express h^(k, rj) in terms of n^(k). With this, Eq. l|Hnjl becomes 
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where again x = krj and xq = krjQ, and we have evaluated the time integral using approximation (|B9|I . Here we have 
also neglected a term of the order of (€^/xq) j£+3(a;o), which in principle is of the same order in the above expression, 
but is always subdominant once we perform the integral over k. Since the power spectra for the E polarization are 
parity even, only the parity even part of the 11='= auto-correlator (Eq. H37() ') contributes to the expression for Cf^ 
derivable from Eq. Ht)3|) . Again we present here only the effect coming from the helical part of the magnetic field, 
using Eq. (|A9|I we find [xq is renamed x) 
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The corresponding equation for C^gy can be found in Ref. jSj- There, a somewhat different approximation than ours 
has been used for the time integral. 

For UA > —2, the integral over x is dominated by the upper cutoff, xd — korjo. Using the approximation ljB10|l . 
we obtain 
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The result for Cj^^^ is obtained upon replacing ha by ns and flA by Qs (more precisely the factor ( ""^^^ ) has to 



be replaced by T'^{^^'^^) and the factor (n^ - l)/(n^ + 4) by ns/(ns + 3)). For -3 < ua < -2, using HB7|I . we obtain 
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Again the E polarization power spectrum from the symmetric part of the magnetic field spectrum is obtained upon 
replacement of ua by ns and flA by fls- Similar evaluations like the ones presented in the previous paragraph show 
that 
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2. B type polarization 



Like for E polarization, the integral solutions for B polarization in the case of tensor perturbations are given in 
. In the approximation of instant recombination we have 
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where 



f,{x) + 2 



jt{x) 



2 



X 



for ^ > 1 . 



With Eq. H60() we can write the above integral in terms of the tensor sources 11^ (k) 
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where we have again used approximation ljB9|l . Like for the E polarization, in this case also it is the parity even part 
of the magnetic source, /(fc), which contributes to the Ci. Eq. (|63(l takes the form 
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Note that within our approximation, for £ » 1, C'^^JL^ ^{A)i ■ "^^^^ is also the case for C^J^^ and C^J^^ , see @. 
Evaluating the integral using expressions (|B10|) and (|B7|I . for the different ranges of the spectral index ua, we obtain 
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Again, the contributions from the symmetric part are obtained by replacing ^Ia by fls and ua by ng, up to factors 
of order unity and we find 
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Within our approximation, which is better than a factor of 2, we have 
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From ordinary inflationary 



jnCf' for gravity waves, which is comparable to our findings. 



Temperature and E polarization cross correlation 



The symmetric part of the source term, /(fc), can only induce parity even CMB correlators. Besides the power 
spectra for temperature anisotropies and E and B type polarizations analyzed in the previous subsections, it can also 
source the cross-correlation between temperature anisotropy and E polarization. In order to evaluate this contribution, 
we have to substitute into Eq. (|63ll the integral solutions for the tensor mode Eqs. (|68|l and H82|) . to obtain: 
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We can evaluate this integral using (|B7() . and we find, 
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(93) 



for -3<nA< -3/2. 



(94) 



(95) 



In this case also, the contribution from the symmetric part of the magnetic field spectrum to the Q-E correlator is 
always larger than this helical part. 



VI. CMB CORRELATORS CAUSED BY MAGNETIC FIELD HELICITY 



If the source (or the initial conditions) have no helical component, (n+(k)n+(k')) = (n^(k)n~(k')), the above 
correlators are the only non-vanishing ones. However, as soon as the tensor magnetic source spectrum has a helical 
contribution (see Eq. (jSHl) 

(27r)3g(fc) EE -^(n+(k)n+*(k) - n-(k)n-*(k)) ^ o , 

the parity odd CMB power spectra are non zero. This has been observed first in Q, where the vector contributions 
have been calculated. Here we compute the gravity wave contributions. We need again to evaluate Eq. Taking 
into account that the gravity waves components h^(k) are directly proportional to the source components (Eq. H60|) ). 
and considering the parity of the radial functions (Eqs. (|67l 1811 IS5|l 'l 



it is clear that cross correlations between temperature and B polarization Cf^, and between E and B polarization 
Cf^^ cannot vanish, since they are given by momentum integrals of g{k). Using the expression of the tensor integral 
solutions 9^ , Ef (jH^J and Bf (jHU , we can calculate the power spectra Cf^ and Cf^ . 



A. Temperature and B polarization cross correlation 



For temperature and B polarization cross correlation we obtain after integrating over time 

e'^ r"dfcfc2:^^±^M^^±iM(n+(fc)n+*(fc)-H-(fc)n-*(fc)) . (97) 

Jo X2 

The antisymmetric source function g{k) is given in Eq. H49|l . and the integral over k can be calculated using (|B7(I . 
Note that g{k) depends on both the spectral indices nA and rig, and we will have to evaluate the integral dividing 
the two cases + ^ — 2. We finally arrive at 
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(98) 



Independently on the spectral indices, t'^Cf^ is always negative for positive Aq. 

In this case of temperature and B polarization cross correlation, we have computed the spectrum (|98|l also numeri- 
cally, in order to test the reliability of our analytical estimation. The amplitude of the numerical result is bigger than 
the analytic one by a factor of two or less, so within the error we estimated for our approximations (see Appendix B). 
We expect this to be one of the worst approximations due to the relatively slow convergence of J dxJi+3{x)Je+i{x)/y^. 



B. E and B polarization cross correlation 



Following the same procedure as in the previous paragraph, we can evaluate the E and B polarization cross 
correlation created by the helical part of the magnetic field. Using the formula 1)63(1 . we get: 
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In the case ua + ns > —2, the mtcgral in x — krjQ is divergent, and we need to evaluate it using approximation l)B12|l . 
which gives: 
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It is not possible to assign a precise value to the variable xd — Vq^d, because of the unavoidable incertitude in the 
estimation of the magnetic field damping scale, which depends on the amplitude of the magnetic field and is therefore 
smeared out over a certain range of scales. Therefore, we expect that the presence of the term sin(2xD) most probably 
leads to a considerable suppression in the amplitude of the E — B cross correlation term. 

For nA + ns < —2, the momentum integral in Eq. (|99|l is dominated by the second term in the square brackets, 
and in order to perform the integration, we need to distinguish two different cases: For —4 < + ng < —2, the 
exponent of x is still positive, so that we have to use the approximation given in Eq. ljB12|l . A further distinction is 
therefore necessary, since the dominant term in approximation (|B12|) depends on whether the exponent is above or 
below 1 as discussed in the Appendix. 
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Both contributions are suppressed by the presence of the two terms sin(2£^) and sin(2a;D) since, usually one averages 
over band powers in £ (for the second case) and also xd is not a very sharp cutoff but has a certain width, as mentioned 
above (for the first case). 

If —6 < nA + ns < —4, the second term in the integrand of Eq. H99(l still dominates, but since the exponent of x is 
now negative, the integral converges and we can make use of approximation ljB7(l . 
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This result is not suppressed by oscillations. 



VII. DISCUSSION AND CONCLUSIONS 



In this paper we have computed CMB anisotropics due to gravity waves induced by a primordial magnetic field. 
We have mainly concentrated on the effects of a possible helical component of the field. Magnetic fields induce scalar, 
vector and tensor perturbations which are typically of the same order. In this sense the tensor contribution can be 
regarded as an order of magnitude estimate for the full contribution. 

As it has already been found in Refs. the C^'s are proportional to 

The first term is (jf-)'^ — lO^^*' (S/lO^^Gauss)^, hence for a primordial magnetic field of the order of i? ~ 10^^ to 
10~® Gauss we would expect to detect its effects in the CMB anisotropy and polarization spectrum. Here B — Bk,^ — 
B\{Xk]j)"^^ is the maximum value of the i?-field which is always the field at the upper cutoff scale l/fc^i which we 
also denote by Bk^ ■ 



In Eq. (|105|l fls stands for fls or flA and in the above expressfon for , n stands for ua or ns depending 
on which contribution we are considering. The second term represents the fogarithmic build up of gravity waves, 
In^ (zi-n/zcq) ~ 660 to 3100. Here the first value corresponds to magnetic field generation at the electroweak phase 
transition, Tin — 200 GeV and the second value represents a possible inflationary generation at Tin — 10^^ GeV. For 
scale invariant spectra, ua = ns — —3, the right hand side of Eq. (|105|1 gives roughly the amplitude of the induced 
CMB perturbations. 

Taking into account the pre-factor 2*-y/7r/9, scale invariant magnetic fields produced at some GUT scale, T ~ 10^^ 
GeV have to be of the order of i? ~ ,B ~ 10"-'^° Gauss to contribute a signal on the level of about 1% to the CMB 
temperature anisotropics and polarization. 

If the initial magnetic field is not scale invariant, the scales ku and rjo suppress the results by factors of l/(fc£)?7o) 
and £/{kDr]o) which are much smaller than unity. Note that the reference scale A introduced in Eqs. does not 

enter in the final results at all, since it is of course arbitrary. 

As already discussed, the damping scale kjj is given by kj^^ ~ VAl^iTcq) ~ x 0.35 Mpc, and va is the Alfven 
velocity, v\ — (i3)^/(47r(p + p)) for the magnetic field averaged over a scale larger than the damping scale. Clearly, 
VA ^ 10^"^ so that B does not induce density perturbations larger than 10^^. Therefore, the damping scale is 
of the order of 1 kpc or less. The latter value is reached for maximal magnetic fields which are of the order of 
(B) ~ 10~^Gauss. On the other hand ao(ryo — ??dec) — Vo is simply the angular diameter distance to the last scattering 
surface, which has been very accurately measured with the WMAP satellite r]o = cIa = 13.7 ± 0.5 Gpc. So that 
koVo 10'' or even larger, depending on the magnetic field amplitude. 

Our results differ somewhat, but not in a very significant way from the results obtained in Ref. jgi]. Since our 
magnetic field spectra are either scale invariant or blue, the induced spectra £'^Cg are also either scale-invariant or 
blue. They grow towards large £. It is therefore an advantage to choose £ as large as possible. However, in our 
calculations we have not taken into account the decay of gravity waves which enter the horizon before decoupling. 
Our results are therfore correct only for £ < r^o/jydcc ~ 60. To be on the safe side, we choose ^ = 50 in our graphics. 

In Fig. n we show ^^C(^-)P at £ = 50 for the different quantities (temperature anisotropy, E and B polarization 
and correlators) as a function of ua with ns fixed to 2 and —2.99. We show the absolute value of the correlator in 
units of 

( ^kVin^ ( 5}A ^ 10-10 ( V 

\nj \z,qj ~ Vio-^Gaussy ' 

and 

^A^s ,2 f fin\ _ 10 f y f Bk„ y 

172 y^^^J l^l0-9Gaussy VlO^^Gaussy ' 

Note that the correlators c'"^'^^ and C'^^^ are always negative and have to be subtracted from cj^^^^ which is of the 
same order of magnitude or larger since fls > ^a and ns < ua- For the limiting case, ~ flA and ns — nA, the 
presence of an helical component in the magnetic field spectrum can in principle cancel the effect of the symmetric 
part on the CMB. In that very particular case, the signature of the presence of a magnetic field will appear only 
through the parity odd correlators. 

From Fig. ^ it is clear that only for n^.s < —2 and ^Ia — ^s ^ 10^^, the effect on the CMB will be of the 
order of a percent or more. In Ref. |0| it has been shown that for ns > —2, magnetic fields with Bx > lO^^^Gauss 
over-produce gravity waves on small scales which is incompatible with the nucleosynthesis bound, for A ~ 1 Mpc. 
Here we require < 10^^ Gauss so that flB remains a small fraction of the radiation density throughout. Then 
Bx = iJfej-, (Afcu)"^""'"'^-' <C Bk„ for n > —2. Therefore, by keeping B^j^ sufficiently small, we automatically satisfy the 
bound derived in Ref. The result is most interesting for the window of —3 < < nA ^ —2 and ^Ia — ^ 10^^, 
which requires Bk^ — B^^ ^ 10^^"Gauss. Especially, if magnetic field helicity is causally produced which implies 
ns = 2 and nA = 3, this effect cannot be observed in the CMB since the parity violating terms are suppressed by 
about 15 orders of magnitude (see lines in the lower right corner of the bottom panel of Fig. 

In Fig. 121 we show the ratio C^(^A)fJ ^{A)t "^-S ~ function of n^. Again, we are mainly interested in the 

part of the graph with —3 < nA < —2, where this ratio raises from the order unity to about 10^. Hence if a close to 
maximal helical magnetic field, with a spectrum not too far from scale invariant, — 3 < ns < nA < — 2 is produced in 
the early universe, it is more promising to search for its parity violating terms than for the parity even contributions. 

We can conclude that helical magnetic fields with a spectrum close to the scale invariant value, —3 < ns — nA ^ —2 
and close to maximal amplitudes on small scales, B^^ > 10~^° Gauss can lead to observable parity violating terms 
C^^ and C^^ in the CMB. Such magnetic fields might in principle be produced during some infiationary epoch where 
the photon is not minimally coupled or via its coupling to the dilaton (see |24l |25| for various proposal of magnetic 
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FIG. 1: On the top panel we show the amphtudes of the parity even correlators, ^'^Cj^^^ (solid, black), 

i'^Cf^^y^ (dotted, red) and ^'^C^^^^' (dashed, blue) as a function of the spectral index ha for £ = 50. 

The logarithm of the absolute value of i^C'^^^'' is shown in units of (QA/^r)'^ ln^(zin/zcq). We do not 

plot ^^C^'^j^' which equals f-^C\^^f^ within our approximation. The spikes at ua = —2 for l^C\^^^ 
and at ua ~ —3/2 are not real. They are artefacts due to the break-down of our approximations at 
these values. 

On the bottom panel we show the corresponding parity odd correlators, ^^C|^^^' (solid, black), ^C\'^^ 

(dashed, red) in units of (SlA^^s/^^r) ln^(ziii/zcq) for ris = —2.99 and ns = 2. In this last case, only 
the allowed range nA > ns = 2 is plotted. Again the spike at ua = 1 for ns = —2.99 and the 
precipitous drop at ua = — 1 in ^^Cj^j^', are due to the limitation of our approximation close to the 
transition indices. 



field production during an inflationary phase). However, so far no concrete proposal has led to ns,A — ^3, nor to the 
creation of a helical term. As we have shown, the effect is largely suppressed and clearly unobservable for causally 
produced magnetic fields, e.g. , during the electroweak phase transition or even later. 

Nevertheless, our calculation also demonstrates the effect of parity violating processes during inflation which may 
lead to a non- vanishing helical component of gravity waves, 7i 7^ 0, see Eq. (I57|l . In this case the above calculation 
can be trivially repeated and will result in non-vanishing parity violating CMB correlators, C^^ 7^ and C^^ 7^ 0. 
We think that already this remark, together with our knowledge that at least at low energies, nature does violate 
parity, should be sufflcicnt motivation to derive experimental limits on these correlators. 



FIG. 2: We show the ratio of the correlators, 
ns ~ —3 as functions of the spectral index ua 
shown in units of Qa/^s < 1- The spikes visible 
due to our relatively crude approximations. 



"a 

Cf^Vcf^^ (solid, black), and C^^^^'/Cr^' for 
for ^ = 50. The logarithm of the absolute value is 
at certain values of the spectral index ua are mainly 
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APPENDIX A: THE SOURCE FOR GRAVITY WAVES 



In this appendix we present some details on how to compute the gravity waves source functions /(fc) and g{k). The 
first step is to evaluate the two point correlator of the magnetic field stress-energy tensor H2I|I : using Wick's theorem 
(122(1 and definition after a longish but simple calculation we obtain 

(ry(k)TL(k')) - i^5(k-k') J d^piSip)Si\k-p\Wu-miSjn.-i^U^)m) 
+ ((5im ~PiPrn){Sjl - (k - p) j (k - p)/ )] 

-A{p)A{\'k - p|)[ej;fej„,.pt(k - p),. + eimfejigPfCk ~ p)g] 

+iS{p)A{\k - p\)[ejmriSil - PiPl){'k - p)r + ejlg{Sim -piPm)(k - p)g] 

+iAip)Si\k-p\)[eM{Sjm - (k^)j.(k^)„)pt + e™/(<5,7 - (k^)j.(k^);)p/] } 

+ ■■■S^-j +---Sim . (Al) 

The isotropic tensor spectrum in the case of a magnetic field spectrum without hclicity term is derived in ^ . Here 
we concentrate on the source terms which contain the helical part of the magnetic field spectrum. 

By acting with tensor projector on I|A1|I . we find expressions H33|) and 1341) for the symmetric and helical parts of 
the source spectrum. Taking into account that the angle /3 = k • (k — p) = . fc~P7 ^ .^^g ^g^^^ rewrite the two 

expressions which contain A{k) in the form 

pA/,,^ _ 8 f ,3^A/„\ Anu „l^ ^ ' - pj) 



r{k) = — ^ / d'pA{p)Ai\k-p\) ,'\,'' (A2) 



^k^ - 2fcp7 + p2 "V k^-2kpj + p- 

The contribution to f{k) from S alone is computed in Ref. 4]. There one finds 



(A3) 



^ - (4^ J - + ^ ) + P - 2kp;+p^ ) • (^4) 

We can now substitute the power law Ansatz (|6I7() for 5* and A in these expressions and try to calculate the integrals. 
The integration over 7 = k • p is elementary, using 

1 Q + 2 

d-f (k^ + p^ - 2kp^)~ = -— —{k'^+p^-2kp^)~ 

kp[a + 2) 

d77" ik^+p' - 2kpj)^ = - /" J fc^+P^ - 2A:p7)^ + , , / ^77""' (fc' - 2A;p7)^(A5) 

kp[a + 2) kp[a + 2) J 

This last integration by parts has to be performed in the worst cases three times, reducing the power m of 7 from 3 
down to 0. 

Since we are integrating 7 over the interval [—1, 1], we get a series of m + 1 terms of the form 

(fc+p)«+2n±|fc_p|a+2n 

{kpT ' ^ ' 

with n = 1, 2, ...(m + 1). To evaluate the integral over p, we can expand those terms using the binomial decomposition 
(1 + a;)" = 1 + ax + a{a — l)x'^ + ■ • •• Since, in general, the value of the exponent a is not an integer, we need 
to truncate the series somewhere, which is well justified only if x ^ 1. To achieve this, we split the integral into 
two contributions, = Jq + J^'^ ■ In the first term p/k < 1, while in the second k/p < 1, which allows us to 
approximate Eq. (|A6|) truncating the binomial series at the second term, 

(k+p) ~\k-p\ - |2ap"-ifc+ia(a-l)(a-2K-3fc3 p>k ^^^> 



and 



(fc+p)" + ifc-pr 



2p°' + a{a-l)p"-^P p>k 



(A8) 



We then perform the integration over p. For each contribution we keep only the terms which, depending on the value 
of the spectral index, may dominate the result. So, we finally obtain, for k < kjj 
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where the coefficients are given by the magnetic field amplitudes at scale A: 
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Aa 



2{2ns + 3) 

2A3 
3(2nA + 3) 
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Bl 



Bl 

2r(ii^) 



3 [ns + nA + 2) 



Bl 



Bl 



2Y{!1A±±^ 



nA - 1 
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nA-i 
ns + 3 



(Afc)"«+"^+2 



(A9) 
(AlO) 
(All) 

(A12) 



(A13) 
(A14) 
(A15) 



The first part of fik), which is the contribution from the symmetric part of the magnetic field power spectrum, has 
been taken from The singularities at n^, = — 3/2 respectively and at ns + = — 2 are removable. 



APPENDIX B: USEFUL MATHEMATICAL RELATIONS 



Integrals of Bessel functions 



In Section [3 we use approximate solutions for the three integrals 



doc ^ (^0 -x)'^ 



(xo - x) 



dx —^jeXxo - X) 



(Bl) 



These integrals are solvable only by numerical method. However, the aim of this paper is to give an approximate 
analytic result. In this appendix we therefore derive and test analytic approximations to the above integrals. To 
achieve this, we first modify them slightly, in order to make them solvable analytically. Then, we adjust the result 
obtained in this way by comparing it with the exact numerical integration. 

Let us concentrate, as an example, on the first integral. We first perform a variable transform to y — xq ~ x. The 
integration boundaries then become and xq — Xdoc- Below, we derive an approximation for 



j2ixQ ~ y) jeiy) 



xo-y 



dy 



Since Bessel functions change on a scale Ay ^ 1, this approximation is good for the integrals in Eq. (jBl|) if x^cc < 1- 
After the integration over x in Ea. ljBl|l we have to perform an integration over k. For £ fixed, this integral is either 



dominated by the contribution art krjQ = xq = £ ot a,t the upper cutoff, fc/j. For the integrals which are dominated 
at Xq = krjo ~ i, the inequaUty Xdcc < 1 is equivalent to ^ ~ — 60a;doc ^ 60. In some cases, however, our integral 
over k is dominated at the upper cutoff kjj with Tydccfcu 3> 60 and of course also Tyofc^ S> 60. Since for i ~ 60, the 
dominant contribution to the integral comes from y 60, our inaccuracy of the boundary will not invalidate the 
approximation also for this case. 

The approximation in the upper boundary of the integral, xq ~ Xdcc — xq makes us miss the characteristic decay of 
fluctuations on angular scales corresponding to £ > 60. 

To make the first integral in Eq. (jBljl solvable analytically, we now modify the powers of y and xq — y- Taking into 
account that the spherical Bessel Function j,y{x) has its maximum value at x v, we make the attempt: 

"\^ hix)je{xo-x) ^ TT j-^" ^^ J5^2{x) Ji+^{xo~x) 







X {xo — xY 2 7q x3/2 (2:0 — x)^/ 



TT 



2_ n J5/2{xa-y) Je+i {y) 
2V5W0 xo-y y 



5 V 5£ 



dx ^' _l . (B2) 



(B3) 



For the last equality, we have used 6.581.2 of [2€ 



b-u N-i T / ^ X / ^ 2'' ^ (-l)"r(6 + p + m)r(6 + m) ,^ ^ 

dxx [a-x) Jp{x)Jq(a - X) = — > — ] nn [b + p + q + 2m) Jb+p+q+2ni{a) , 

aq ^-^ m\ I (b) I [p + m + 1) 

(Re {b + p) > 0,Req> 0) (B4) 

and the recurrence relation J^-i{x) + J^+i{x) = — (9.1.27 of T^), keeping only the highest order terms in £. We 
can now compare this approximated analytic result with an exact numerical integration. Since the analytic result is 
again a Bessel function divided by a power law, it has a maximum at xq ~ and its envelope has a power law decay 
for Xq > £. This two characteristics are very well reproduced by the numerical result, which however decays somewhat 
faster; it turns out that a better approximation is 



Xq 



j2ix) je{xo - x) 1 /3£ Je+3ixo) 







X (xo - a;)2 3 V 2 -3 y ' 



To estimate the goodness of our approximation, let us now take into account the integration over /c, as in Eq. 163|l . 
What we are finally interested in is (Eq. ifTTll l 



CLXq Xq g 



, / \ 271 A+3' 

nA - 1 / Xq 
HA + 4: \Xd 



(B6) 



As already discussed in the main text, this integral is always convergent and dominated by the contribution around 
Xo — £ : we should therefore make sure that our approximation is good around that value. We have that for £ — 30, 
our approximation underestimates the numerical result by about factor of two; for £ = 40, the error reduces to 15%, 
and is always smaller for larger values of £. 

Fig. El shows the numerical result for the integral in l|B5(l (green, dotted line), together with its analytical approxi- 
mation (the right hand side of Eq. ljB5|) , blue and long dashed) and a numerical evaluation of the same integral when 
x^cc is not set to zero (red, solid). For small values of £ (in the left hand panel of Fig.O £ = 50), Eq. ljB5|l is a good 
approximation in the region xq ~ £. However, if ^ > 60 setting x^cc — > causes a large overestimation of the result. 
In the right hand panel of Fig. El it is shown that, for £ — 100, the difference between the integral with lower bound 
and the one with lower bound Xdcc is of more than a factor of ten. Consequently, as already stated before, we can 
rely on all our approximations only for £ < 60. 

We proceed now to evaluate integral ljB6p . Since xd = korjo > 10^, for £ < 60, integral (|B6|I can be calculated in 
the limit xd ^ 00, using formula 6.574.2 of ['26j |: 

r(6)r(2±2^') 

dx Jp{x)Jq{x)x-'' = -, r ) ^ (B7) 

^ nbj- f -P+q+b+l \ p f p+q+b+l \ p / p-g+b+O 

(Re(p + q + 1) > Re6 > 0) . 








FIG. 3: In both panels, as a function of xo'- the green dotted line shows the numerical value of the integral 
in <B5ll . the blue, long dashed line shows the analytic approximation (right hand side of Eq. IBSt 'l. and 
the red, solid line shows the numerical value of integral IBSII if Xdcc is not put to zero. All these functions 
are squared, and multiplied by Xq: this gives us an indication of the result, after the integration over xo, as 
stated in Eq. IIB6II . In the left panel £ — 50, in the right panel £ = 200. First of all, we note that it appears 
clearly that the value of the integrals is dominated at xq — £, and that the function goes to zero quicker 
than Xq'^, which justifies our approximation xd ^ 00 and the use of formula IB 71 Secondly, we note that 
for ^ = 50 and xq ~ £, our approximation (blue, long-dashed) is good for both the integrals. However, if 
£ = 200, the approximation overestimate the correct numerical result by about a factor of ten. 



This approximation is used for example in Eqs. (|74l I75|) . 

With the same procedure we can approximate the second integral of Eq. ljBl|l . for which we find {£ < 60) 



^^JMM^^IM^JI^. (B8) 
dec ^ (^0 ^ X) 3 V 2 Xq 

This approximation underestimates the numerical result with an error of about 40% for i = 30, which reduces to 20% 
at ^ = 60. In this case also, the integral over xq is convergent, and we can proceed as before. 

The situation is different for the third integral of Eq. (jBl|l . In this case, the numerical result is approximated by 
the following function {£ < 60): 

"° , J2{x) . , ,1 [2 J,+^{xo) 

dx ]i(xo-x)~-\ . (B9) 

xdoc ^ 3 V 5 ^ 



It is clear that if we insert this function in an integral like (jB6|l we cannot perform the limit xd ^ 00 since this 
integral is dominated at the upper cutoff. Consequently, we need a good approximation for the behavior of the 
integral for large values of xq — > xjj. In this case, we no longer require our approximation to be accurate at xq ~ £, 
but we concentrate on its behavior for high values of xq, which will dominate in the integral over xq. Fig. ^ shows 
the approximation for i = 30, which overestimate the numerical result by an error within 1%. 

We also have to evaluate the integral over xf, dxo of the square of ljB9|l . which we encounter in two different cases. 
The first (see Section IV B(l is of the kind J^" dxx^jf^x). For p < this integral converges and we may evaluate it in 
the limit xd — > 00, in which it is of the form ljB4|) . For p > and xd 3> £'^, the integral can be approximated using 
the asymptotic expansion of j£{x) for large arguments Je(x) ^ -\/2/(7r.T) cos[.t — {2£ + l)7r/4]. Approximating 
the oscillations by a factor of 1/2, we obtain 



dxxPjfix)- dxxPjf(x) ~ i P^' , (BIO) 
Ji2 \ iln^ , p = 



1000 



100 




0.01 



FIG. 4: We plot the value of integral 1IB9II squared and multiplied by Xq as function of xo, for £ — 30. The 
green, dotted line represents again the numerical result (a^dcc ~^ 0), and the blue, long dashed line is the 
analytic approximation. In this case the slope is positive, and hence the integral dxo/xo of this function is 
dominated by the upper cutoff. 



For the second case, dx xPJe{x)j£-^^i{x), which we encounter in Section IVII we use again the large argument 
approximation for the Bessel functions, for x ^ 



Je(x)Ji+i(x) ~ —cos fx - (2^+1)^) cos fx - (2£ + 3)v) = — cos f x - (2^+1)^) sin fx - (2^+1)^ 
TTX V 4/V 4/7rxV 4/V 4 

(_l).+i 



TTX 

— sin ( 2x 

TTX 



■ cos(2x) 



(Bll) 



so that for p > 



dx x^ Ji{x)Ji^i[x) ~ ' 



dxxf-i cos (2x) ~ -^^^ (x^^^ sm{2xD) - ^^""^ sm{2f)^ 



(B12) 



In the limits to which we have restricted ourselves, we always have xd 3> (."^ ■ Consequently, the dominant contribution 
in the last expression can be given either by the first term in the bracket, if p > 1, or by the second term, if p < 1. 
Numerical checks show that the approximation is good for p > 1, but it is rather poor in the second case, p < 1. Since 
we shall not be very much interested in this case, we do not go any further in this work. 
When evaluating expression (|B7|I . we often also use 



r(2x) = r(x)rfx + - 



r(x + a) 
r(x + 6) 

(see Eqs. (6.1.18) and (6.1.47) of 0). 



0{x' 



,a-6-l\ 



for X 1 



(B13) 



2. Recurrent Relations for spherical Bessel Functions 



We use several recurrence relations for spherical Bessel functions in our derivations, most notably 



(B14) 



and 



-jiix) -j'iix) = ji+iix) . 



(B15) 
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